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We study the effect of gluon number fluctuations (Pomeron loops) on deep inelastic scattering 
(DIS) in the fixed coupling case. We find that the description of the DIS data is improved once gluon 
number fluctuations are included. Also the values of the parameters, like the saturation exponent 
and the diffussion coefficient, turn out reasonable and agree with values obtained from numerical 
simulations of toy models which take into account fluctuations. This outcome seems to indicate the 
evidence of geometric scaling violations, and a possible implication of gluon number fluctuations, in 
the DIS data. However, we cannot exclude the possibility that the scaling violations may also come 
from the diffusion part of the solution to the BK-equation, instead of gluon number fluctuations. 
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The mean-field dynamics of the high-energy dipole-proton scattering is described by the BK-equation 1]. Phc- 
nomcnological ansatze for the dipole-proton scattering amplitude T(r, x) (where r is the transverse dipole size and x 
the Bjorken-variable) inspired by the BK-equation have led to quite successful descriptions of the HERA data. The 
T-matrix following from the BK-equation shows within a restricted kinematical window, which increases with collision 
energy, the geometric scaling behaviour | 11, i, T(r,x) = T(r 2 Q 2 s (x)), where Q s (x) is the saturation scale, which 
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Q. seems well supported by the HERA data |f| . The correction to the solution outside the restricted window, the "BK- 



diffusion term" , violates the geometric scaling 0, S 13 and depends on the variable ln(l/r 2 Q 2 (x))/\/IJq^Y. Iancu, 
Itakura and Munier (IIM) [fj have shown that the "BK-diffusion term", giving a substantial amount of geometric 
scaling violations, is needed in order to accurately describe the experimental HERA data. The exponent A of the 
i—i, saturation scale, Q 2 s {x) ~ (xq/x) , is known at NLO [7(, A ~ 0.3, and agrees with the values extracted from fits to 
HERA data. 

Recently, there has been a tremendous theoretical progress in understanding the high-energy QCD evolution beyond 
• the mean field approximation, i.e. beyond the BK-equation. It has been understood how to include discreteness and 
fluctuations of gluon numbers (Pomeron loops) in small- a; evolution [1, 0, [13, El, E3 ■ After including these elements, 
t— ( ■ the evolution becomes stochastic and one has to distinguish between the event-by-event amplitude T(r, x), which 
corresponds to an individual gluon number realization, and the physical amplitude (T(r,Y)}, which one obtains by 
^ — . averaging over all individual realizations [9(. At very high energy, the discreteness effect decreases the exponent A as 
compared to BK- value and the gluon number fluctuations, i.e., the averaging over all events to calculate the physical 
t^"- , amplitude, replaces the geometric scaling resulting from the BK-equation (in the "wave front" region) by a new 
scaling the diffusive scaling, namely, (T(r,Y)} is a function of a single variable ln(l/r 2 Q 2 s {x))/ V 'DY ', where D 

is the diffusion coefficient. The value of D determines the rapidity above which gluon number fluctuations become 
important, Y > Yd = 1/D, which is the case when the fluctuation of the saturation scales of the individual events 
becomes large, in formulas, when the dispersion er 2 = 2((p 2 (Y)) — (p s (Y)) 2 ) = D7» 1, where p s (Y) = \n(Q 2 (Y)/Q 2 ). 
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At high energy, such that a 2 1, it has been shown that fluctuations do strongly modify measurable quantities |13l.ll4l|. 
(A more detailed presentation of the recent theoretical pro gres s is given in R efs.flal while the most recent studies on 
Pomeron loops based on toy models can be found in Refs. JLsl . ll7L [HI, H, H3, l2ll.l22l [23|J3] .) 

In this work we elaborate, in a quite approximative way, whether the HERA data [25| do indicate any possible 
implication of gluon number fluctuations. The coupling is kept fixed throughout this work. We proceed in the 
following way: We use for the event-by-event amplitude several models, the GBW model [26|, the IIM model ((J and 
a model which is close to the theoretical findings for T at very large energy (see Eq. ([7])). For the averaging over 
all events we use the high-energy QCD/statistical physics correspondence [3], i.e., a Gaussian for the distribution of 
p s (Y) — \n(Ql(Y) / Qq) . Moreover, assuming that the DIS cross section shows diffusive scaling in the HERA energy 
range, we have used the "quality factor" method of Ref. [27] to get an estimation for the value of A, in a model- 
independent way. The procedure we use in this work is always based on approximations and, therefore, can at best 
give hints on a possible implication of gluon number fluctuations in the HERA data. 
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After including fluctuations in the way described above, we obtain from the analysis of the HERA data values 
for the exponent A and the diffusion coefficient D which are quite independent of the ansatz for the event-by-event 
amplitude. Also the model-independent approach gives a similar value for A. We find A ~ 0.2 which is smaller than 
the value from the BK-inspired models (no fluctuations), A ~ 0.3, and the decrease is in agreement with theoretical 
expectations. For the diffusion coefficient we find a sizeable value, D ~ 0.35. Surprisingly, this value is very close to 
the values found for D in numerical simulation of the (1 + 1) dimensional model and of evolution equations in 
QCD [28] (approximations to Pomeron loop equations O EL [HI ) m the nxe d coupling case. The sizeable value of D 
may indicate a possible involvement of fluctuations in the HERA data since Y > Yd = 1/-D for rapidities at HERA. 

We observe that after including fluctuations the description of the HERA data is improved for all models we have 
used for the event-by-event amplitude. In the case of the GBW model, which exhibits pure geometric scaling, after the 
inclusion of fluctuations, which lead to a violation of geometric scaling, a much better description is obtained, namely, 
X 2 /d.o.f = 1.74 without and x 2 /d.o.f = 1.14 with fluctuations. The situation seems to be similar with all event-by- 
event amplitudes which show geometric scaling. In the case of the IIM model, which contains already the geometric 
scaling violating BK-diffusion term, the inclusion of fluctuations also improves, however less than in the GBW case, 
the description of the HERA data; x 2 /d.o.f = 0.983 before and x 2 /d.o.f = 0.807 after including fluctuations. 1 The 
outcomes seem to tell us that violations of geometric scaling are required for an accurate description of the HERA data. 
The improvement of the description of the HERA data together with the very reasonable values for the parameters 
discussed above seem to indicate that gluon number fluctuations may be the reason for geometric scaling violations 
in the HERA data. However, we wish to emphasize here that the BK-diffusion term gives similar geometric scaling 
violations as fluctuations and may as well be the reason for the geometric scaling violations in the HERA data. 

This work is organized as follows: In Sec. [II] we show the results for the T-matrix for dipole-proton scattering and 
for the energy dependence of the saturation scale which are obtained in the mean field approximations, i.e., from the 
BK-equation. The results for the same quantities beyond the mean field approximation, or the effects of discreteness 
and fluctuations in gluon numbers on these quantities, are summarized in Sec. IIII1 Finally, we give numerical results 
and discuss a possible implication of the physics beyond the mean field approximation in the HERA data. 



II. MEAN FIELD APPROXIMATION 



In the mean field approximation, the Y-dependence of the T-matrix for a dipole of transverse size r scattering off 
a proton is given by the BK-equation. In the fixed coupling case, the solution to the BK-equation in the saturation 
region, where T ~ 1, is [29] 



T(r,Y) = 1-Co exp[-d(p- Ps (Y)) 2 ] for p - p s (Y) <C 1 , 
while for the front of the T-matrix, where T <C 1 (but not too small) , one finds [2J, |j] 



(1) 



T{r,Y) =C 2 [p- Ps (Y)+C 3 ] exp 



-X s (p-p s (Y)) 



2a X "(X s )Y 



for l«p-p s (y)«2 X "(A s )a s y , (2) 



where have used a s — a s N c /ir, p = \n(l/r 2 Qy) and p s {Y) — \u(Q 2 (Y)/Qq) with Q S (Y) the saturation scale. In 
above equations, the constants Co, C*2, C3 are of 0(1), C\ — —Cp(l — Xo)/N c 2x{X s ) (Cf is the casimir factor in the 
funcdamental respresentation), A s = 0.6275, and xW = 2^(1) — 4>(X) — — A) is the eigenvalue of the BFK1 kernel. 



For the rapidity dependence of the saturation scale, which_segarates the saturated (r l/Q s (Y)) from the dilute 

BK equatii 

Ql(Y) = Ql 



h sepa 

(r <C 1/Q S (Y)) regime, one obtains from the BK equation [l,|j] 

exp[a X '(A s )r] 



(3) 



Note that within the even more restricted window, p — p s (Y) -^i yj^x' (A s )a s y, where the diffusion term in the 
exponent in Eq.Q can be neglected, the T-matrix shows the geometric scaling behaviour, i.e., it depends only on the 
difference p — p s (Y) intead of depending on r and Y separately. At very large r, so that p — p s (Y) 3> 2x"(X s )a s Y, 
the T-matrix exhibits color transparency, i.e., it shows a faster decrease with p as compared to Eq.([2]); T ~ exp[— p\. 



1 The x 2 is defined such that the smallest x 2 gives the best description to the HERA data. 
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lancu, Itakura and Munier Q have used the following ansatz for the T-matrix, 

1 - exp [-a In 2 (brQ s (a:))] , rQ s (x)>2 

2(A i + '°' 2 a"» ) 



7V 1 ' , rQ s (x)<2 



which obviously includes the features of the solution to the BK equation, to compare the theory in the mean field 
approximation with the DIS data. They have used for the saturation momentum the leading Y-dependence of Eq((3} , 
Q s (x) = (xq/x) x , however, with A and Xo being fixed by fitting the DIS data. The constant n = x"{^s)/x'(^s) ~ 9.9 
is a LO result coming from the BK-equation, Nq is a constant around 0.5 and a and b are determined by matching 
the two pieces in Eq.Q at rQ s = 2. 

The "BK-diffusion term" in the IIM-ansatz flU), 



ln(2/rQ s (x)) 



: A Y 

= exp 



ln 2 (4/r 2 Q 2 s (x)) 



2k\Y 



(5) 



which is the quadratic term in the exponent of Eq.Q, does explictly violate the geometric scaling behaviour. We wish 
to emphasize here that, as also shown in this violation seems required in order to get an accurate description of 
the DIS data. Without it, even allowing A s to be an additional fitting paramter, one can not get a better description 
of the DIS data. For further details on the importance of the diffusion term see Ref. @ . 

In this work, we wish to elaborate whether the violation of the geometric scaling may come from gluon number 
fluctuations (Pomeron loops) and not from the BK-equation. As we will see in the next sections, the fluctuations 
do indeed give a similar violation of the geometric scaling and also lead to a better description of the DIS data as 
compared to the case where the T-matrix shows a geometric scaling behaviour. 



III. BEYOND THE MEAN FIELD APPROXIMATION 



To go beyond the mean field approximation one has to include the effect of discreteness and fluctuations of gluon 
numbers [a 0|. After including fluctuations one has to distinguish between the even- by-event amplitude and the 
averaged (physical) amplitude. They can be explained by considering the evolution of a proton from y = up to y = Y 
which is probed by a dipole of size r, giving the amplitude T(r, Y). The evolution of the proton is stochastic and leads 
to random gluon number realizations inside the proton at Y, corresponding to different events in an experiment. The 
physical amplitude, T(r,Y), is then given by averaging over all possible gluon number realizations/events, T(r,Y) = 
(T(r, Y)), where T(r, Y) is the amplitude for the dipole r scattering off a particular realization of the evolved proton 
at Y. In the following we discuss the event-by-event amplitude T(r,Y) and the averaged amplitude T(r,Y). 



A. Event-by-event scattering amplitude 



In a single scattering process, the mean field approximation breaks down when the occupancy of gluons inside the 
evolved proton is low so that the discreteness of the gluon number needs to be taken into account; the number of 
gluons cannot be non-zero and less than one since it has to be discrete. When including the discreteness effect, as 
compared to the results from the BK-equation, the energy dependence of the saturation momentum changes to [8, 9] 



Ql(Y) 



)l exp 



d sX '(X s )Y 1 



2(Ap) 2 x (A s ) 



(6) 



and the piecewise, approximate, shape of the T-matrix at fixed coupling and very high energy reads [8, 9] 



I for p-p s (F)<0 

T(r,Y) = I Ni [p- p s (Y)}e x ^P-P° (Y K for < p - p s (Y) < Ap 
, N 2 e^-^)] for p - p s (Y) > Ap 



(7) 
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where Ni and N2 are irrelevant constants and the front width is Ap ~ (1/A S ) ln(l/a 2 ). The front width cannot be 
larger than Ap which is the distance when the amplitude decreases from its maximal value T w 1 down to the value 
T = 0(a 2 ) where the discreteness of gluon numbers becomes important. The width is formed via diffusion, p—p s (Y) oc 
yj aY, and it requires the rapidity Yf ~ (Ap) 2 /(2x"(A s )a) until it is completed. The event-by-event amplitude given 
in Eq.©, which is formed at Y > Yp, shows, approximately, geometric scaling: T(r, Y) w T(p — p s (Y)). 

The main differences as compared to the mean-field results are: The exponent of the saturation scale in the event- 
by-event amplitude, cf. Eq.© and Eq. ©, is decreased due to the discreteness of gluon numbers. Further the width 
of the front of the event-by-event amplitude is fixed, Ap, instead of increasing with rapidity as in Eq.©. 



B. Physical scattering amplitude 



Based on the relation between high-energy QCD evolution and reaction-diffusion processes in statistical physics 0] , 
the fluctuations in gluon numbers are taken into account by averaging over all event-by-event amplitudes, 

(8) 



(9) 



(T((p-p s (Y)))= J d Ps T(p- p s (Y)) P(p s (Y) - (p s (Y))) , 
where the distribution of p 8 (Y) is, to a very good approximation, a Gaussian [3l| : 



: exp 



ip. - (Ps)Y 



The expectation value of the front position, (p s (Y)), increases with rapidity as (p s (Y)) = ln(Q 2 (Y)/QQ) at high 
energy [9], with Q S (Y) given in Eq. (©. The dispersion of the front at high energy increases linearly with rapidity, 



a 2 = 2[(p 2 s )-(p s ) 2 ]=DY 



(10) 



where D is the diffusion coefficient, whose value is known only for a — > (asymptotic energy) [g, l33 |. Since the 
values of D and the exponent A of the saturation scale, Q 2 {x) = 1 GeV 2 (xq/x) x , see Eq. ((6|), are not known for finite 
energies, e.g. at HERA energy, in what follows we will treat them as free parameters. 

At very high energy, such that a 2 3> 1, the dispersion of the fronts due to the gluon number fluctuations from event 
to event has large consequences on (T(r,Y)): the geometric scaling of the single events T(p — p s (Y)), cf. Eq.©, is 
replaced by a new form of scaling, known as diffusive scaling, namely, (T(r, Y)) is a function of (p — (p s (Y)))/\/DY , 



(T(r,Y))=T(r,Y)=T 



P-(Ps(Y)) 



(11) 



The diffusive scaling is expected to set in at Y > Yd = l/D, which follows from the requirement a 2 ^> 1. 

The goal of this work is to study whether the diffusive scaling behaviour of the dipole-proton scattering amplitude 
in Eq. (|11[) . which is caused by gluon number fluctuations, may be present in the HERA data. As we will see in the 
next section, the fluctuations do improve the description of the HERA data, indicating that the violation of geometric 
scaling seems important for an accurate description of the data. We will discuss whether the violation prefered by the 
DIS data is due to the gluon number fluctuations, which lead to the diffusive scaling (p — p s (Y))/ V DY , or due to the 
BK diffusion term, cf . Eq. ([5]) , which corrects the geometric scaling in a similar way, namely, via (p— p s ( Y) ) / y/2 ax" Y. 



IV. NUMERICAL RESULTS 



Our fit includes the ZEUS data for the F2 structure function, 

F 2 (x,Q 2 )= %- (<j T (x,Q 2 )+a L (x,Q 2 )) , a T ,h{x,Q 2 ) = f dz d 2 r \^ T x(^r,Q 2 )\ 2 a dip (x,r) (12) 

47T Ct em J 

in the kinematical range x < 10~ 2 and 0.045 GeV 2 < Q 2 < 50 GeV 2 (see also [6| for more discussions on the range). 
The upper limit on Q 2 has been chosen large enough to include a large amount of "perturbative" data points, but 
low enough in order to justify the use of the BFKL dynamics, rather than DGLAP evolution. We use in our fit the 
same photon wave functions ipr,L as m Ref.[26(, which are computable in QED, and three quarks with equal mass, 
m q = 140 MeV. We have considered only the ZEUS data because there is a mismatch between the HI and ZEUS with 



model / parameters 


x 2 


X 2 /d-o.f 


x (xlO" 4 ) 


A 


R(im) 


D 




266.22 


1.74 


4.11 


0.285 


0.594 







173.39 


1.14 


0.0546 


0.225 


0.712 


0.397 



TABLE I: GBW model: The parameters of the event-by-event (2 line) and of the physical (3 line) amplitude. 



regard to the data normalization and since only ZEUS has data also in the low Q 2 region, i.e., in the saturation region. 
To fix the parameters we minimize \ 2 = Si( m °del(?,pi, ■■■,Pn) — -F2(i)) 2 /(error(i)) 2 , where the sum goes over the 
data points, p\,...,p n denote the parameters to be found, ^(i) the experimental results for the F2 structure function, 
and for the error of F2, i.e., (error(i)) 2 , we use the systematic error squared plus the statistical error squared. 

The interesting ingredient for us in Eq. (fT2|) is the dipole-proton cross section, ddip = 2ttR 2 (T(r,x)), with 2nR 2 
being the outcome of the integration over the impact parameter. We will use different ansatze for the event-by-event 
amplitude, T(r, x), and the physical amplitude, (T(r,x)), is obtained according to the rules outlined in section fill Bl 
(We wish to note that the ansatze for T(r,x), which are derived/motivated based on perturbative QCD, are used 
to describe also the low virtuality data, Q 2 < 1 GeV 2 , in the fit to the HERA data. In this region non-perturbative 
physics [33I ] is involved which is only approximately given by our ansatze.) In Odip we will use the event- by- event 
amplitude and the physical amplitude in order to study the effects of gluon number fluctuations. In the case of T(r, x) 
there are three free parameters which will be fixed by fitting the HERA data: R ( "radius of the proton" ) and xq and 
A coming via the saturation momentum Q 2 {x) = 1 GeV 2 (xq/x) x . In the case of the averaged (physical) amplitude, 
(T(r,x)}, there is another free parameter, the diffusion coefficient D. 

• Golec-Biernat, Wiisthoff (GBW) model [26|: 



The GBW model 

■ 2 Q 2 s (x) 



T uaw {r,x) = 1-exp 



(13) 



is one of the most simple models which shows geometric scaling, T(r,x) — T(r 2 Q 2 (x)), and leads to a quite 
successful description of the HERA data, as can be seen from Figs. [TJ [5] and the x 2 (error) in Table U (denoted 
by GBW). It is nice to see that the value of the saturation exponent, A ~ 0.285, which is found by fitting the 
HERA data with the GBW model, comes out close to the theoretical NLO results for A 0] • 

Now, using the GBW model as an event-by-event amplitude, we include the effect of gluon number fluctuations 
by averaging over all events via Eq. ([5]). The resulting (T GBW (r, x)), which breaks the geometric scaling, leads 
to a relatively much better description of the HERA data, as can be seen from the comparision of the % 2 values 
and the two lines in Figs. [TJ [21 The large improvement after including fluctuations seems to indicate that 
violations of geometric scaling, and probably even gluon number fluctuations, are implicated in the HERA data. 

It is important to note that the values of the fitting parameters come out reasonable also after including the 
gluon number fluctuations. The value of A becomes smaller after including fluctuations which is in agreement 
with theoretical expectations, as can be seen from the comparison of Eq. ([3|) with Eq. ([6|). Furthermore, the 
value of the diffusion coefficient D is sizeable, and is surprisingly close to the values which have been found 
numerically by solving the (1 + 1) dimensional toy model [l9[ and the a ppr oximate QCD evolution equations (28j 
(they represent an approximation of the Pomeron loop equations (Tol. Illl [T2I]) in the fixed coupling case. Note 
also that the radius of the proton, R, increases somewhat and xq becomes smaller, meaning that Q s < 1 GeV up 
to x ~ 10 -6 , due to fluctuations. Also the reasonable values of the parameters, especially the sizeable value of 
D yielding Yd — l/D — 2.5, in addition to the better description of the HERA data after including fluctuations, 
seem to be in favour of an implication of gluon number fluctuations in the HERA data. 

Iancu, Itakura, Munier (IIM) model [6i], other models and a model-independent approach: 



The IIM model given in Eq. ^ includes the BK-diffusion term, ln(4/r 2 QD/v2k\Y, which explicitly violates 
the geometric scaling. It has been shown in Q that this violation does noticeably improve the description of 
the HERA data in comparison with the GBW model, as can be seen from the much smaller \ 2 value in the IIM 
case in Table |TT] (we always use N = 0.5 in the IIM model). In Ref. @ has been further shown that without 
the BK-diffusion term, although allowing for an additional free parameter A s (one parameter more than in the 
GBW model), the x 2 /d.o.f value does not improve and is close to the GBW value. 



model / parameters 


x 2 


X 2 /d-o.f 


xo (xKT 4 ) 


A 


ii(fm) 


D 


T IIM 


150.45 


0.983 


0.5379 


0.252 


0.709 







122.62 


0.807 


0.0095 


0.198 


0.812 


0.325 



TABLE II: IIM model: The parameters of the event-by-event (2 line) and of the physical (3 line) amplitude. 



Note that the GBW model only after including gluon number fluctuations gives a x 2 /d.o.f value which is 
comparable with the IIM one. This may mean that the violation of the geometric scaling is favoured by the 
HERA data. The violation may come from the gluon number fluctuations or from the BK-diffusion term. 



In the case of the IIM model, 
amplitude 



after including fluctuations, we can give an analytic expression for the physical 



(T IIM (r,y)) 



No 
2a 



In- 



Exp 



crErfc 



ln 2 (6 2 r 2 Q 2 ) 



1 + Erf j 



x Erfc 



4 ^ <r 2 ul Vr 2 Q 2 I 



+ h Hwb*) 



2k\Y 



x Exp 



2kXY 



2k AVer 2 4 + As 



V 



2k\Y 



which can be used in phcnomcnological applications, where Erfc(x) is the complementary error function. Also 
in the IIM case fluctuations do improve the description of the HERA data, however not much, as can be seen 
from the comparable x 2 / d -°- f values for T IIM and (T IIM ) in Table HU This is so because the IIM model does 
already contain the geometric scaling violations via the BK-diffusion term, ln(4 / V 2 Q 2 S ) / V 2kXY, and describes 
accuratelly the HERA data, before including fluctuations. However, note that the diffusion coefficients in case 
of fluctuations and the BK-diffusion term are quite different, namely, D = 0.325 and 2 k A ~ 3.9, respectivelly. 

After including fluctuations, the parameters in the GBW and the IIM case are close to each other. Apart from 
the fact that similar values for D are found in numerical simulations of evolution equations [HI, [28| and the 
decrease of A due to fluctuations is theoretically expected, at least at high energy, the parameters A and D also 
seem to be quite model-independent. Indeed, similar values for A and D would come out also if one uses a model 2 
as suggested by the thoeretical findings at high energy as given in Eq. l[7]). for reasonable values of the proton 
radius, R ~ 0.7 — 0.8 fm. Moreover, the similar value of A coming out of the different models is also supported 
by the following model-independent approach: In case fluctuation are important in the range of HERA data, 



'DY is a function of ln(l/r 2 Q 2 )^ 



Defining a 



one finds the diffusive scaling behaviour 30], i.e., a 1 
"quality factor" C(A, Xq,D) as done in [27], which test the quality of this diffusive scaling in HERA data, we 
find by minimazing it A = 0.215, at least for the input-values 0.01 < D < 0.7 which we have investigated. 

The seemingly model-independent values of the parameters A and D, their agreement with numerically found 
values, and the improvement of the description of the HERA data in all models after including fluctuations, 
seem to tell us that gluon number fluctuations are relevent in the range of HERA data. However, since in 
the case of the IIM model the fluctuations do not improve much the description of the HERA data, one may 
conclude that the BK-equation alone should describe the HERA data and that fluctuations are neglegible in the 
energy range of the HERA data. The intention of this work is to illustrate the possibility that fluctuations may 
be present in the HERA data. 

The situation would become more clear at even higher collision energies as compared to the HERA energy. 
With growing Y, according to the BK-equation the window for the geometric scaling behaviour would increase, 
and the scaling violating term would become less important. On the other side, the small- a; dynamics including 



2 Such a model would be for instance the IIM model with the diffusion variable ln(4/r 2 Q 2 )/\/2 k A Y replaced by ln(4/r 2 Q 2 )(l — X s ) / Ap, 
such that the new model interpolates between the three regions of Eq. (JTJ) and shows the geometric scaling behaviour. The constant Ap 
is given by Eq.j7]l. We use in Ap a small value for a, a = 1/15, which is the value required such that the exponent of Q 2 in Eq.j3]l agrees 
with experimental or NLO results, A ~ 0.3. With this input, we find for R = 0.8 fm, the following results: A = 0.235 and D = 0.58. 
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gluon number fluctuations leads to a more clear diffusive scaling behaviour with increasing Y. The forthcoming 
LHC may tell us more whether geometric or diffusive scaling is more appropriate for the description of the 
observables in the LHC energy range. 

Throughout this work the coupling is kept fixed. As mentioned above, the (f + 1) dimensional model in [lj|, 
which accomodates Pomeron loops, gives similar values for D as our analysis for a fixed coupling. However, it 
has been recently shown that if allowing the coupling to run in this toy model pol ] then the effects of gluon 
number fluctuations can be neglected up to energies far beyond the HERA and LHC energies. We plan to 
extend our work by the running coupling in order to see whether the HERA data can tell something about the 
running coupling and whether the prediction of the toy model remains valid also in the QCD case. 
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FIG. 1: The Fi structure function versus x at different values of Q 2 . The experimental points are the latest published data 
from the HI (squares) and ZEUS (dots) collaborations (25|. The solid line represents the result of the averaged GBW fit and 
the dashed line the result of the GBW fit to the ZEUS data for x < 1CT 2 and 0.045 GeV 2 < Q 2 < 50GeV 2 . The data points 
at lowest Q 2 values, 0.045, 0.065 and 0.085 GeV 2 , are not shown here although they are included in the fits. 




FIG. 2: The same as in Fig. [T] but for larger values of Q 2 . Note that we show in this figure our results up the highest Q 2 
although our fit is performed inluding only the data for Q 2 < 50 GeV 2 . 



